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1.  SUMMARY 

The  research  completed  under  the  contract  F44620-75-C-0064 
"Research  in  Seismology"  during  the  period  1  Jemuary  1978  - 
30  June  1978  falls  within  the  two  broad  topics  of  (1)  Seismic 
source  mech2uiisms  and  (2)  Seismic  wave  propagation.  The 
specific  research  tasks  within  each  category  have  broad 
applicability  to  the  problem  of  discrimination  of  earthquakes 
from  iinderground  nucleeir  explosions,  with  particular  emphasis 
on  the  Asian  continent. 

Source  mechanism  investigations  included  continued  detailed 
studies  of  earthquakes  in  eastern  Turkey  from  the  stemdpoint 
of  field  and  teleseismic  observations  and  of  precursory 
phenomena  CNabelek  and  ToksOz,  1978a, b)  and  a  new  look  at 
the  earthquake  mechcinism  auid  stress  drop  from  the  standpoint 
of  a  model  fault  with  barriers  (Aki,  1978) .  Wave  propagation 
work  included  studies  of  attenuation  mechanisms  in  rock  (Johnston 
and  Cheng,  1978} ,  attenuation  and  shear  wave  delay  beneath  the 
Mid~Atlantic  Ridge  (Muller  et  al. ,  1978),  emd  the  further 
application  of  a  technique  for  simultaneous  inversion  of 
surface  wave  attenuation  €uid  phase  velocity  information  to  data 
from  several  continental  emd  oceanic  paths  (Lee  cund  Solomon, 

1979) . 

Details  of  these  studies  are  given  in  the  preprint  and 
abstracts  in  the  following  sections.  A  cumulative  list  of 
publications  completed  under  the  contract  is  also  included. 
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2.  SEISMIC  SOURCE  MECHANISMS 


SOURCE  PROPERTIES  OF  THE  1976  EARTHQUAKE  IN  E.  TURKEY 

NABELEK,  J.L.  AND  TOKSOZ,  M.N.,  DepartmenC  of  Earth  and  Planetary 
Sciences,  Massachusetts  Institute  of  Technology,  Cambridge, 

HA  02139 

The  Eastern  Turkey  earthquake  on  24  November  1976  (Ms  *  7.3)  provided 
good  data  for  testing  the  validity  of  source  models  by  comparing 
geological  observations  at  the  source  and  seismic  observations  at 
teleselsmlc  distances.  A  steeply  dipping  fault  was  traced  for  55  km 
with  an  average  strike  of  N70*V.  The  motion  was  essentially  right-> 
handed  strike  slip  with  an  average  displacement  of  2.5-3  m.  Water 
level  marks  along  the  shore  of  Lake  Van  provided  an  estimate  of  the 
regional  tilt. 

The  fault  plane  solution  obtained  from  P-wave  polarities  and  the  surface 
wave  radiation  pattern  Is  well  constrained  with  a  strike:  N73*U,  dip: 
78*S  and  slip  angle:  4*.  The  width  of  the  fault  was  about  25  km, 
based  upon  static  dislocation  theory  and  shape  of  the  surface  wave 
spectra.  Seismic  moment  calculated  on  the  basis  of  static  and  dynamic 
source  models  was  1.4  x  10^^  and  0.7  x  10^^  dyne-cm  respectively. 

Slight  disagreement  between  two  values  may  be  due  to  frequency 
dependence  of  the  moment,  aselsmlc  creep  or  aftershocks.  The  stress 
drop  was  35  bars. 


Presented  at;  Seismological  Society  of  America  Annual  Meetina 
J^ril,  1978 
Sparks,  Nevada 
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Eastern  Turkey  Earthquake  of  24  Novenber  1976  - 
A  Comparative  Study 


Dlst 


J.  Nabelek  and  H.N.  ToksBz 


The  24  November  1976  earthquake  In  the  Van  Province  of  eastern 
Turkey  was  a  large  (Mg  ~  7.3)  strike-slip  event.  A  comparison  of  this 
earthquake  with  the  1976  Guatemala  event  can  provide  a  better  under¬ 
standing  of  the  dynamics  of  strike-slip  faulting.  Geological 
observations  in  the  source  area  of  the  Turkish  earthquake  are  well 
documented.  Surface  expressions  of  the  fault  could  be  traced  clearly 
for  55  km.  Average  displacements  along  the  dislocation  vectors  were 
3  -  2.5  m.  The  dip  was  almost  vertical.  Water  level  marks  along  the 
shore  of  Lake  Van  provided  an  estimate  of  the  regional  tilt.  A  few 
physical  precursory  phenomena  were  also  observed.  Seismic  moment  and 
source  mechanism  were  determined  both  from  field  observations  In  the 
source  region  and  seismic  records  at  teleseismlc  distances.  The 
agreement  Is  excellent  for  the  orientation  of  the  fault.  However, 
seismic  moment  determined  from  surface  wave  data  (7  x  10^^  dyne-cm)  Is 
about  2  times  smaller  than  the  static  moment  estimate.  The  reverse 
was  observed  for  the  Guatemalan  earthquake.  The  geometry  of  the  two 
events  Is  different.  The  Turkish  fault  zone  Is  relatively  short  and 
wide  (55  km  and  25  km  respectively),  with  large  displacements  (2.5  m). 
Vaulting  In  the  Guatemalan  event  was  about  250  km  long  with  average 
surface  displacements  of  about  1.0  m.  Differences  between  the  source 
parameters  of  the  two  earthquakes  may  be  due  to  differences  in  stress 
regime,  age  of  the  fault,  or  thickness  of  the  brittle  layer. 


Presented  at;  International  Symposium  on  the  February  4,  1976 
Guatemalan  Earthquake  and  the  Reconstruction 
Process 
May,  1978 
Guatemala  City 


STRESS  DROP  AND  ENERGY  RELEASE  IN  AN  EARTHQUAKE 
ESTIMATED  BY  THE  USE  OF  THE  BARRIER  MODEL 

.  Kellcl  Akl,  (Dept,  of  Earth  &  Planetary  Sciences, 

Mass.  Inst,  of  Techn. ,  Cambridge,  MA  02139} 

So  far,  the  stress  drop  In  an  earthquake  has 
been  estimated  by  assuming  a  fault  plane  without 
any  obstacles,  mismatches,  or  barriers.  On  the 
other  hand,  observations  of  slip  along  a  fault  at 
the  surface,  observations  of  short'period  seismic 
waves  near  a  fault  and  at  faffield,  and  scaling 
relations  for  small  earthquakes  in  various  seis¬ 
mic  regions  suggest  that  the  fault  plane  with 
distributed  barriers  which  remain  unbroken  after 
rupture  propagation  may  be  more  realistic  than  a 
smoothly  ruptured  fault  plane.  By  numerical 
experiments  on  spontaneous  shear-crack  propaga¬ 
tion,  Das  and  Akl  (1977)  demonstrated  the  mech¬ 
anical  feasibility  of  such  a  model  of  a  fault 
plane  with  distributed  unbroken  barriers. 
Application  of  the  barrier  model  to  actual  data 
leads  to  a  disturbing  conclusion  that  the  stress 
drop  in  major  earthquakes  might  have  been  under¬ 
estimated  by  an  order  of  magnitude.  The  model 
also  requires  the  apparent  stress  to  be  an  order 
of  magnitude  greater  than  earlier  estimates. 

Then,  we  face  another  disturbing  result  in  that 
the  Gutenberg-Richter’s  magnitude-energy  formula 
may  be  underestimating  seismic  energy  by  an 
order  of  magnitude.  The  estimated  high  ambient 
tectonic  stress  is  more  consistent  with  the  fric¬ 
tional  strength  of  rocks  measured  in  laboratories. 
Such  a  high  stress,  however,  cannot  be  maintained 
by  plate  boundary  forces  such  as  slab-pull  and 
rldge-push,  but  requires  a  basal-shear  due  to 
convection  in  the  upper  mantle  as  discussed  by 
Hanks  (1977). 


Presented  at:  American  Geophysical  Union 
Spring  Annual  Meeting 
April,  1978 
Miami  Beach,  Florida 


3.  SEISMIC  WAVE  PROPAGATION 


THE  ROLE  OF  FLUID  FLOW  IN  ATTENUATION 
David  H.  Johnston 

C.U.  Cheng  (both  at:  Dept,  of  Earth  and 

Planetary  Sciences,  Massachusetts  Institute 

of  Technology,  Cambridge,  1L4  02139) 

The  relative  importance  of  Biot  type  flow  and 
"squirting"  mechanisms  In  determining  the  overall 
attenuation  in  saturated  sandstones  at  various 
frequencies  and  pressures  is  studied.  Elastic 
properties  for  a  multiphase  medium  are  found 
using  the  technique  described  in  ToksBz  et  al. 
(1976).  These  moduli  may  be  complex,  introducing 
attenuation.  Biot  or  inertial  flow  predicts  an 
f2  dependence  for  the  attenuation  coefficient 
at  frequencies  below  that  of  Poisseuille  flow 
breakdown  and  an  f^'^  dependence  above.  The 
contribution  due  to  this  mechanism  is  found 
using  a  formulation  of  Biot  theory  adapted  by 
Stoll  and  Bryan  (1970).  "Squirting"  flow 
involves  the  relative  change  in  volume  between 
connected  cracks.  The  attenuation  is  found  by 
introducing  a  complex  fluid  bulk  modulus  deter¬ 
mined  by  the  relaxation  time  for  the  flow.  The 
critical  frequency  for  this  mechanism  is  in  the 
sonic  range.  Another  related  mechanism,  viscous 
shear  relaxation,  is  of  importance  only  at  very 
high  frequencies  (>1  MHz)  for  pore  fluids  of  1 
cp  viscosity. 

Theoretical  calculations  of  each  mechanism 
for  the  Berea  sandstone  are  carried  out  based 
on  ultrasonic  attenuation  measurements  as 
functions  of  differencial  pressure.  It  is 
concluded  chat  for  this  case,  fluid  flow  plays 
a  secondary  role  relative  to  ocher  attenuation 
mechanisms  such  as  friction  across  cracks  and 
grain  boundaries.  Ac  these  ultrasonic 
frequencies  and  relatively  high  pressures  prior 
to  pore  collapse,  Biot  flow  plays  a  more 
important  role. 

Stoll  and  Bryan,  197Q,  J.  Acoust.  Soc.  Am..  47. 

1440-1447. 

ToksiJz  et  al.,  1976,  Geophysics.  41.  621-645. 


Presented  at:  American  Geophysical  Union 
Spring  Annual  Meeting 
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ATTENUATION  AND  TRAVEL  TIME  DELAYS  OF 
SHEAR  V7AVES  BENEATH  THE  MID-ATLANTIC 
RIDGE  FROM  TRANSFORM  EARTHQUAKES 

James  L.  Muller 
Jennifer  Lynn  Hall 

Sean  C.  Solomon  (all  at:  Dept,  of  Earth 
and  Planetary  Sci.,  Massachusetts 
Institute  of  Technology,  Cambridge, 

MA  02139) 

For  earthquakes  on  ridge-ridge  trans¬ 
form  faults,  the  azimuthal  variation  of 
body-wave  differential  attenuation  and 
travel-time  delay  may  be  used  to  explore 
the  shallow  structure  of  the  transform 
and  of  the  intersecting  ridge  crest  seg¬ 
ments.  For  the  13  February  1967  earth¬ 
quake  on  the  Gibbs  transform,  SH  waves 
show  much  higher  attenuation  for  paths 
beneath  the  ridge  crest  to  the  northwest 
than  for  other  directions  [Solomon,  1973] 
cuid  the  S  wave  travel-time  delays  also 
show  a  strong  azimuthal  dependence  [Du- 
schenes  and  Solomon,  1977].  We  have  re¬ 
peated  these  two  sets  of  measurements  for 
a  number  of  additional  earthquakes  on 
mid-Atlantic  ridge  transforms,  including 
the  Gibbs,  Oceanographer,  Romanche,  and 
St.  Paul's  transforms.  Within  the  uncer¬ 
tainties  in  the  respective  data  sets,  no 
clear  direction  dependence  is  apparent 
for  azimuths  near  the  transform  strike 
for  either  SH  differential  attenuation  or 
S  wave  delay  for  any  of  these  additional 
events.  Thus  velocity  and  Q  variations 
beneath  young  seafloor  in  general  are 
either  not  as  pronounced  or  not  as  exten¬ 
sive  as  those  inferred  from  the  results 
of  the  13  February  1967  event,  and  either 
that  segment  of  the  Gibbs  transform  or 
that  particular  earthquake  source  may  be 
anomalous. 


Presented  at:  American  Geophysical  Union 
Spring  Annual  Meeting 
April,  1978 
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ABSTRACT 

^A  formalism  for  the  simultaneous  inversion  of  surface 
wave  phase  velocity  and  attenuation,  previously  developed 
for  Love  waves,  is  extended  in  this  paper  to  Rayleigh  waves. 

The  simultaneous  inversion  technique  permits  the  specification 
of  the  intrinsic  dispersion-attenuation  relation  that  arises 
from  linearity  and  causality,  and  takes  full  account  of  the 
dependence  of  surface  wave  phase  velocity  and  on  the 
and  imaginary  parts  of  an  anelastic  earth  structure.  The 
formalism,  including  resolution  analysis  and  extremal  inversion, 
is  applied  to  combined  Love  and  Rayleigh  wave  data  sets  for  a 
tectonically  active  and  a  stable  continental  path  and  to 
Rayleigh  wave  data  for  a  stable  oceanic  path.  The  depth 
to  the  low  Q  zone  is  60^0  km  for  the  central  Pacific, 

80^20  km  for  western  North  America,  and  130+30  km  for  east- 
central  North  America,  q"^  within  the  low-Q,  low-velocity  zone, 
however,  is  greater  beneath  western  North  America  than  beneath 
the  central  Pacific;  the  low-Q  zone  beneath  east-central  North 
America  need  not  be  a  low  velocity  zone  at  frequencies  above 
1  Hz.  'Kie  surface  wave  data  cannot  distinguish  among  several 
possible  intrinsic  dispersion-attenuation  relations  for  the 
upper  mantle.^ 
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INTRODUCTION 

There  are  several  advemtages  to  the  simultaneous  treatment 
of  phase  velocity  emd  Q~^  data  in  the  inverse  problem  for 
surface  wave  dispersion  and  attenuation.  The  two  data  sets 
are  coupled  because  both  depend  on  the  complete  auielastic 
structure  of  the  earth.  More  fundamentally,  if  linearity  holds 
for  seismic  waves,  a  reasonadale  assumption  at  the  small  seismic 
str2d.n  amplitudes  (less  than  10*^).,  then  there  is  a  predictable 
relation  between  local  body  wave  phase  velocity  and  attenuation. 
In  particular,  the  principle  of  causality  leads  to  an  anelastic 
dispersion  in  solids  which  depends  in  detad.!  on  the  magnitude 
and  frequency  dependence  of  attenuation  (Lomnitz,  1957; 
Putterman,  1962;  Jeffreys,  1965;  1975;  Carpenter  and  Davies, 
1966;  Randall,  1976;  Liu  et  al.,  1976;  Kanamori  and  Anderson, 
1977) .  Eventtial  isolation  of  the  physical  mechanisms  which 
control  attenuation  and  anelastic  dispersion  in  the  mantle  cem 
result  only  by  slmulitaneovis  consideration  of  phase  velocity 
and  data.  In  am  earlier  paper  (Lee  amd  Solomon,  1978; 
hereinafter  called  Paper  I) ,  we  presented  a  formalism  for 
simultaneous  inversion  of  phase  velocity  and  Q  ^  for  Love 
waves,  including  the  effect  of  anelastic  dispersion,  amd  we 
applied  the  formalism  to  data  from  western  North  America. 

An  essentially  equivalent  treatment  for  Love  waves  without 
anelastic  dispersion  has  been  derived  independently  by  Silva 
(1978) .  In  this  paper. 
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we  extend  the  inversion  formalism  of  Paper  I  to  Rayleigh 
waves,  auid  we  apply  it  to  combined  Rayleigh  and  Love  wave 
data  sets  in  western  and  eastern  North  America  and  to 
Rayleigh  wave  data  in  the  central  Pacific. 

An  approximate  correction  to  surface  wave  phase 
velocities  or  to  normal  mode  periods  to  account  for 
euielastic  dispersion  has  been  given  by  Carpenter  emd 
Davies  (1966),  Randall  (19760 »  and  Liu  et  al.  (1976).  Using 
this  correction  and  the  assumption  of  Q  ^  independent  of 
frequency  over  the  seismic  wave  band,  Anderson  et  al.  (1977) , 
Anderson  and  Hart  (1976)  and  Hart  et  al.  (1976,  1977) 
inverted  corrected  normal  mode  periods  to  derive  frequency- 
dependent  velocity  models  for  the  earth.  An  important 
result  of  their  work  is  the  removal  of  the  'beiseline 
discrepancy'  between  earth  models  derived  from  short-period 
body  wave  data  and  models  derived  from  long  period  free 
oscillations.  The  simultaneous  Inversion  procedure  presented 
in  this  paper  differs  formally  from  the  approximate 
correction  technique  and  leads  to  somewhat  different  results 
even  for  the  velocity  structure  except  in  the  limit  of  a 
narrow  frequency  band  for  both  data  euid  earth  models  of 
interest  (Paper  I;  Lee,  1977). 

In  the  sections  to  follow,  we  first  extend  to 
Rayleigh  waves  the  formalism  of  Paper  I  for  simultaneous 
inversion  of  phase  velocity  and  Q“^  data  to  obtain  a  complex 
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frequency  dependent  earth  model.  The  theoretical  basis  for 
the  formalism  is  the  generalization  of  Haskell's  (1953) 
matrix  treatment  to  2melastic  structures  by  Schwab  and  Knopoff 
C1971,  1972} .  Resolving  length  analysis  emd  extremal 
inversion  (Lee  amd  Solomon,  1975)  are  applied  to  Love  and 
Rayleigh  wave  data  sets  over  stable  and  tectonically  active 
portions  of  North  America  euid  to  Rayleigh  wave  data  in  the 
Pacific  (Mitchell  et  al. ,  1976) .  Of  interest  are  both  the 
differences  in  the  results  between  simultaneous  inversion  and 
previous  separate  inversions  of  velocity  or  Q”^  data  sets  emd 
the  inqplications  of  the  different  anelastic  structural  models 
for  lithosphere  and  asthenosphere-  for  the  various  paxhs 
saspled.  ' 

FORWARD  PROBLEM 

The  problem  of  surface  wave  propagation  through  perfectly 
elastic  multilayered  media  can  l^e  treated  by  the  matrix 

formulation  of  Haskell  (1953) .  In  each  layer,  with  boundary 

\ 

conditions  of  the  free  surface  and  of  continuity  of  stress 
emd  displacement  at  the  interfaces,  a  set  of  eq..ations  hold: 


(  Xj  +  2pj) 


2 

"3  't  *3  -  “j 


(1) 
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where  7J  «  -L  +  -L 
3x“^  32^ 

3*  3ij;.  3^.  in,. 

Ua  «  — i  +  -r^  ,  w .  «  — i - i 

■*  3x  32  ^  32  3x 


where  x  (propagation  direction)  and  2  are  the  hori2ontal  and 
vertical  auces,  respectively,  Xj,  and  ^ure  Leune  parameters 
and  density  in  the  j-th  layer,  ipj  euid  c^j  are 
potentials  of  the  elastic  field  of  the  jth  layer,  and  u j ,  v j , 

Wj  are  the  displacements  in  the  x,  y,  and  2  directions.  For 
am  anelastic  (or  viscoelastic)  medium,  the  wave  equation  (1) 
and  the  solution  have  the  same  form  in  the  frequency  domain  as 
for  an  elastic  medium  except  that  the  elastic  mod\ilus  is 
replaced  with  a  complex  quantity  according  to  the  'correspondence 
principle'  (Christensen,  1971).  The  Fourier  trauisform  of 
equation  (1)  with  complex  modulus  is 


E*(a))  Vj  Fj  »  Pj  0)2  Fj(x,2,o>)  (2) 

where  Fj(x,2,o))  »  /  F.,(x,2,t)  e^“^  dt 

2iT  —  J 

where  E*j(a))  represents  either  +  2Uj  or  ji ^ ,  Fj 


.4 
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represents  ,  i|)^  or  ,  emd  is  the  Fourier  trwsfom  of 
Fj.  (u)  depends  upon  frequency  in  general.  The  frequency 
dependence  can  be  specifically  defined  by  the  creep  function 
(or  relaxation  function)  of  the  medium.  The  solution  to  equation 
(2)  for  Fj  «  Vj  (for  Love  waves)  with  the  Isoundary  conditions 
mentioned  ^d}ove  is  given  in  Paper  I.  In  analogous  fashion, 
Haskell's  matrix  formulation  may  be  extended  to  lossy  media  for 
Rayleigh  waves  as  well  as  Love  waves  by  implementing  complex 
velocities  amd  a  dispersion  relation  between  the  real  auid 
imaginaury  parts  of  the  intrinsic  velocity. 


INVERSE  PROBLEM 

The  phase  velocity  amd  attenuation  of  surface  waves  on 
a  multilayered,  amelastic  eaurth  au:e  obtained  from  the  roots  of 
the  caiq>lex  dispersion-attenuation  functions  (Schwad)  and 


1971)  fj 

(Love) 

J 

and 

^R 

(Rayleigh) : 

"y 

dj) 

»  0 

fRlTf, 

a*. 

dj)  -  0 

i 

-  1 

r  2, 

e  •  • 

,  m; 

j  “  0,l,2,...,n 

(3) 


where  T^,  c^,  auid  c^  are  ith  period  and  Love  and  Rayleigh 
wavo  phase  velocities  and  Uj,  6j  and  dj  are,  respectively, 
the  P-  and  S-wave  velocities,  density  auid  thickness  for  the 
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jth  layer  (j>0  is  the  water  layer  for  Rayleigh  waves  over 

Xi  R 

oceanic  paths) .  The  velocities  c  ,  c  ,  a*  amd  3*  are  complex 
quantities . 

The  Inverse  problem  consists  of  finding  am  amelastlc 
(complex)  earth  model  from  given  observational  data  pairs, 
phase  velocities  and  phase  attenuations.  We  staurt  with  am 
initial  amelastlc  earth  model  amd  a  set  of  observed  dispersion- 
attenuation  data  pairs.  The  phase  velocity  and  attenuation  for 
the  initial  complex  model  are  then  calculated  theoretically 
by  Haskell's  (1953)  method  at  the  period  of  each  observed  data 
pair.  The  generalized  Haskell  formulation  is  for  a  flat  earth, 
whereas  the  observations  are  for  a  spherical  earth.  The  flat- 
to-spherlcal  transformation  of  Biswas  and  Knopoff  (1970) ,  as 
asiended  by  Schwab  and  Knopoff  (1971)  to  Include  anelasticlty, 
is  used  for  sphericity  corrections  for  Love  waves.  For  Rayleigh 
waves,  a  similair  transform  is  given  by  Schwab  and  Knopoff  (1972) , 
though  such  a  treatment  is  difficult  to  use  in  a  computer  code. 

In  this  study.  Bolt  and  Dorman's  (1961)  empirical  correction  has 
been  used  to  account  for  sphericity  emd  gravity  for  con^lex  phase 
velocity.  Although  the  empirical  correction  has  been  tested  only 
for  the  perfectly  elastic  case,  it  is  probably  sufficiently 
accurate  for  our  purposes.  North  2md  Dziewonski  (1976)  ioproved 
this  correction,  but  we  did  not  adopt  their  improvement  because 
only  a  minimal  change  is  expected  at  the  periods  used  here  (<100 
sec) . 

From  the  dispersion-attenuation  calculations  we  also  obtain 
the  partial  derivatives  of  the  complex  phase  velocity  c  with 
respect  to  each  conqplex  parameter  p  of  the  layered  earth 
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model.  Then  the  initial  perturbation  equation  can  be  written 
relating  the  desired  parameter  corrections  to  the  differences 
between  the  corresponding  theoretical  and  observed  phase 
velocity  values: 

AcJ'®  -  -  AP4  (4) 

^  3Pj  ■' 

where  the  repeated  Indices  isply  svumnation  for  n  layers.  A 

similar  equation  for  each  period  ceui  be  formed. 

Because  the  physical  significemce  of  a  coo^lex  quantity 

is  more  easily  understood  by  decoiqposltion  into  real  and 

imaginary  parts,  we  write  equation  (4)  as  two  real  equations, 

0 

rather  than  one  complex  equation: 
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R  L  R 

where  c^'  and  C2'  are  the  real  and  imaginary  parts  of 

T  D 

c  '  and  and  P2  are  the  real  and  imaginaury  peurts  of  p. 

It  is  assuiaed  that  all  the  dissipation  is  due  to 
imperfect  elasticity.  By  requiring  the  density  to  be  real 
we  ignore  the  possibility  of  losses  due  to  in^rfect  Inertia 
(Anderson  2md  Archanbeau,  1964).  For  Love  waves  we  can  write 
equation  (5)  in  the  form: 
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(6) 


where  and  real  and  imaginary  parts  of  the  shear 

wave  velocity  in  the  jth  layer.  For  Rayleigh  waves, 
equation  (5)  can  be  written  as 
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(7) 


Since  one  or  both  of  A8j^  and  A82  depend  upon  frequency  (as  do 
AU]^  and  Aa2  for  Rayleigh  waves),  the  right  hand  sides  of 
equations  (6)  and  (7)  should  be  standardized  at  a  single 
reference  frequency  for  inversion  (see  below) .  Generalizing 
equation  (4)  to  m  complex  observations  and  using  matrix 


notation  gives: 


.L,R  L,R 
A  X, 


(8) 


where  is  an  m  x  1  matrix  of  differences  between  observed 
and  predicted  phase  velocities  and  attenuations/  is  an  m  x  n 
matrix  of  partial  derivatives/  euid  is  eui  n  x  1  matrix  of 
perturbations  to  the  starting  anelastic  eeurth  model.  The 
elements  of  and  x^  are  real  2x1/  2x3  aund  3x1 

'Vj  A/  % 

p  R  R 

matrices/  respectively  and  the  elements  of  y"/  ^  auid  ;|c  are 
real  2x1/  2x5  and  5x1  matrices/  respectively. 

The  partial  derivatives  of  phase  velocity  with  respect  to 
shear  velocity  tre  obtained  by  implicit  function  theory  (Schwab 
and  Knopoff /  1972) : 


agL/R  ^ 

36*  36*  3c^'® 

-  .!!5  /  !!r 

3o*  3a*  3c** 


(9) 


«d:ere  f^/  given  in  Paper  1/  and 


^(o)-(l)j.(2) 


p(n-2)_(n-l)^(n) 

_{n-2)  (n-l)-(n) 

F  F  T 


even  n 
odd  n 


are  the  dispersion-attenuation  functions  for  Love  euid  Rayleigh 
waves,  and  where  T^°^,  T^"^/  F^^^/  and  are' given  in  Schwab 

and  Knopoff  (1972) .  The  partial  derivatives  required  to  evaluate 
equation  (9)  are  given  for  Love  waves  in  Paper  I.  For  Rayleigh 


waves,  define 


(n-1) 


p(l)p(2)|r(3)  _^(n-2)p(n-l) 


Then 


ij.(o)  ^(n-l)^(n) 


even  n 


(o)  (n-l)_(n) 

T  A  T  for  odd  n 


For  even  n. 
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V  p(1)f<2) 


ac® 


^  -  T^°^F<l)p(2)_ 


+  T^°^F^^^F^^^...-a£li\..ir(n-l)  (n) 
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From  equation  (11) ,  only  four  of  the  elements  of  the  2x3 
real  matrix  in  equation  (6)  are  independent  for  Love  waves 
and  only  six  of  the  elements  of  the  2x5  real  matrix  in 
equation  (7)  are  independent  for  Rayleigh  waves.  Thus  the 
matrices  can  be  completely  specified  from  the  real  and 
imaginary  parts  of  3c^'^/36*and  3c®/3cF  and  from  3 c^'^/3 p. 

RESOLUTION 

As  in  Paper  I,  we  follow  the  treatment  of  Der  and 
Landisman  (1972)  for  parameter  resolution  for  a  two-parameter  (5^/82) 
earth  model.  The.  linear  combination  of  data  used  as 

the  estimator  of  a  desired  parameter  in  layer  Jc  is 
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r  'v^'® 
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(12) 


where  is  a  coefficient  to  be  determined,  N  is  the  number 
of  data,  M  is  the  nuoher  of  layers  in  the  model,  and  where 


3x 


L,R  ' 

j 


(13) 


normalized  by  the  layer  thickness  in  kilometers. 

The  three  quantities  to  be  minimized  ave  (Der  and 
Landisman,  1972) 

1)  the  variance  of  the  desired  variadsle 

2)  S]^,  the  resolution  for 

•“L  R 

3)  82 1  the  dependence  on  the  undeslred  varied) le  Xj^' 
for  the  same  layer. 


± 
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If  we  assume  that  the  observational  errors  are 
independent 


var  f  var 


3-k  ^  *3 


(14) 
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Where  dj  is  the  layer  thickness  and 
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These  three  objectives  cam  be  acconplished  by  minimizing  the 

function 

\ 


*  ain5  var  +  cosCsinn  (s^'^)  +  cos£cosn  (s^'^)  ^ 

+  2e(E^*-l)  (16) 

where  6  is  a  Lagrange  multiplier.  The  parameters  C  and  n 
(0  ^  ^,n  £  "’f/2)  are  adjusted  so  eis  to  balance  the  three 
desired  minimizations.  As  C  is  increased,  the  approximation 
to  the  delta  function  becomes  poorer,  the  varlauices  of 
become  smaller  and  separation  between  and 
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Xj^  '  improves . 

Resolution  analysis  should  strictly  be  conducted  with 
density  as  a  third  independent  parameter,  since  density  is 
included  in  the  inversions  presented  in  this  paper.  With  no 
additional  data,  such  as  from  higher  modes,  the 

resolution  for  6^  emd  ^2  would  be  worse  than  for  the  examples 
shown  below. 

DISPERSION-ATTENUATION  RELATIONS 

Dispersion-attenuation  relations  can  be  given  in  two 
different  forms.  In  the  freqtiency  domain,  for  a  given  frequency 
dependence  of  Q,  Kreuners-Krfinig  relations  provide  this  relation. 

A  frequency  dependence  or  independence  of  Q  can  be  specified  by 
a  superposition  of  Individual  relaxation  mechemlsms  or  visco¬ 
elastic  elements,  or  by  empirical  observations.  In  the  time 
domain,  Boltzmann's  after-effect  equation  will  provide  a  complex 
modulus  emd  a  dispersion  relation  for  a  given  creep  f\inction  or 
relaxation  function.  A  creep  (or  relaxation)  function  can  similarly 
be  given  by  superposing  the  creep  (or  releucation)  function  of  each 
viscoelastic  element  using  a  distribution  function,  or  by 
eiiq>irical  observations.  No  matter  which  procedure  is  used, 
the  same  dispersion-attenuation  relations  are  given  for  the 
same  frequency  dependence  of  Q.  We  discuss  the  two  categories  of 
frequency  independent  euid  frequency  dependent  Q  separately. 

Frequency  independent  Q.  Various  attempts  have  been  made  to 
explain  the  apparently  nearly  constant  Q  in  the  seismic  frequency 
band  (Futterman,  1962;  Lomnitz,  1957;  Aziml  et  ^. ,  1968;  Liu 
et  al.,  1976).  Futterman  (1962)  and  Azlmi  et  al.  (1968)  derived 
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dispersion  relations  in  the  frequency  domain: 


X(w) 


/ 


Cj^cii 


XqM 
1  + 


(Futterman) 

(17) 

(Azimi  et  al. ) 


where  x(*o)  is  the  attenuation  coefficient,  or  the  imaginary 

part  of  the  complex  wave  nxmber,  and  C. ,  x  /  xV  const£uts. 

X  o 

Equations  (17)  give  the  following  dispersion  relations: 


V((i»)/Vo, 


2Ci  V.  _i 

- - —  Inoj]  (Futterman) 


,-l 

— — ln(Xi<i»)l  (Azimi  et  al. ) 
ir(l-Xi  “  ' 


(18) 


where  v<»  is  the  non-dispersive  limit  of  v  and  where  x^^  is 
chosen  so  that  X (u)  is  almost  linear  in  some  finite  frequency 
rimge  of  interest  0  ^  ^  “lim  (^i®i  «t  al.  used  the  value 

10~^  sec  for  the  constant  Xj^)  • 

For  attenuation 


Q"^(u)-  J 


2CiV, 


2x  V 
*0  0 


1  +  Xj^w 


Since  Xj^  is  very  small  (10"^), 
nearly  identical. 


(Futterman) 


(Azimi  et  al. ) 


these  two  expressions  are 
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Lomnitz  (1957)  and  Liu  et  al.  (1976)  instead  derived 
dispersion  relations  in  the  time  domain.  The  creep  function 
was  given  as  follows: 


q  ln(l  at)  (Lomnitz) 


^  -t/T 

C  /(1-e  '  )dT/T  (Liu  et  al. ) 

^1 


(19) 


where  q/  a,  and  c  are  constants  and  t  is  a  relaxation  time. 
These  creep  functions  give  the  following  dispersion  relations 


il  +  q/2rY+  In  |j 
1  +  C/2  In  (uT^) 


(Lomnitz) 
(Liu  et  al. ) 


(20) 


where  y  is  a  constwt,  and  1/Q  is  approximately  constant  at 
the  veLlue  qV2  or  CTf/2,  respectively.  From  (18)  euid  (20), 
«rtien  Q“^  2L  Qo^  small. 


v((iij^)/v(u2) 


1  “l 

—  ^ 

’^Qo  2 


This  is  a  good  approximation  for  various  attenuation  laws  of 
constant  or  nearly  constant  Q. 
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Frequency  dependent  Q.  Jeffreys  (1958)  modified  equation  (19) 
of  Lomnitz  (1957)  to  the  Jeffreys-Lomnitz  law,  which  also 
represents  the  empirical  relation  of  Andrade  (1911) .  The 
creep  function  for  the  Jeffreys-Lomnitz  law  is 

V 

♦(tl  -  q/v  [(1  +  at)  -  11  (22) 

where  q  cind  a  are  constants  and  v  is  a  number  between  0  amd  1. 
As  V  goes  to  zero,  (22)  tends  tow2urd  the  original  Lomnitz  Law. 
In  the  frequency  domain,  the  creep  function  (22)  gives  the 
following  dispersion  relations 


Q“^(a>)  »  q  a'*(v-l)l  u>“'^  sin  ^ 

0 

»  II  +  5 —  cot 

2  2 


(23) 


(Jeffreys  and  Crampin,  1960;  Jeffreys,  1965,  1975) .  A  number 
of  suggestions  for  the  proper  value  of  v  have  been  made  by 
Jeffreys  (1/4,  1/5)  and  Andrade  (1/3).  Lamb  (1962),  Strick 
(1967)  and  Azlmi  et  al.  (1968)  considered  the  following 
frequency  dependence  of  the  attenuation  coefficient: 


X(“) 


C2 


1-v 

(I) 


(24) 


where  v  Is  a  number  between  0  to  1 
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The  frequency  domain  approach  gives  the  following 
dispersion  relations  from  (24)  (Paper  I) : 


Q“^(u))  «  2C2 

V(uj)/Va,  »  [1  +  ^  ^ 


(25) 


which  is  identic2Ll  to  (23)  for  suitable  choice  of  C^-  Lamb 
(1962)  used  v  s  1/2  and  Strick  (1967)  and  Azimi  et  al.  (19€8) 
attempted  to  account  for  nearly  constant  Q  with  v  »  0.078 
(Strick)  amd  v  -  0.1  (Azimi  et  al. ) . 

Equations  (23)  and  (25)  give' 
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where  »  Q~^(<*»2^* 

Other  dispersion-attenuation  relations  are  possible  with 
different  frequency  dependence  for  v(w)  and  Q'‘^((d);  for 
instance,  the  relation  for  a  standard  linear  solid  is  given 
in  Zener  (1948) .  Only  relations  (21)  and  (26)  will  be  used, 
however,  in  the  inversions  below. 
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INVERSION  PROCEDURES 

■  For  the  dispersion-atteniiation  relations,  Qg^,  Q,,  »  and  the 
real  body-wave  phase  velocities  S  and  a  are  related  to  ^2' 


and  Qj  such  that 

Q-l  -  26j/6i 


(27  ) 


Q-1  -  202/0^ 

e  -  (ej  +  82^/Si 
a  *  (a^  + 


For  each  of  these  relations,  Sj^,  or  all  B^,  Bj'  ®2 

depend  on  frequency 

6i^C«i)  -  Bij  +  g(»i)^2j 

(28) 

Oij(Ui)  -  o^j  +  g(»i)02j 


a2j  (Wj,) 


Ojj  h(u)j^) 


28. 


where  g  and  h  sure  specified  functions  of  frequency  and 
®lj  ®2j  values  of  8^ ^ and 
at  a  reference  frequency,  respectively.  In  general,  the 
inverse  problem  to  equation  (8)  is  conducted  at  the  reference 
frequency,  the  p£u:tlals  in  (6)  emd  (7)  are  with  respect  to 
^1'  ^2'  ®1'  ®2'  earth  structure  at  any  other 

frequency  follows  from  (28)  (see  Appendix) . 

DATA 

The  data  sets  for  emd  for  North  American  paths 

are  described  emd  tabulated  In  Lee  emd  Solomon  (1975) .  The 

acconpemying  phase  velocity  measurements  are  given  In 

Solomon  (1971).  The  western  North  American  path  crosses 

primarily  the  tectonically  active  Basin  and  Range  physiographic 

province.  The  east-central  North  Amerlceui  path  crosses 

primarily  the  tectonically  quiet  Great  Plains  emd  Central 

-1 

Xiowletnds  physiographic  provinces.  The  Q  data  were  also 
analyzed  from  a  different  viewpoint  by  Solomon  (1972a, b) . 

The  data  set  for  for  the  central  Pacific  Ocean  is 
a  subset  of  that  given  by  Mitchell  et  al.  (1976) .  These 
workers  determined  from  the  records  of  three  earthquakes  at 
WWSSN  stations  distributed  around  the  west  coast  of  America, 
the  Far  East  and  the  Pacific.  The  measurements  represent  a 
weighted  average  of  the  entire  Pacific  Ocean  region.  Such 
an  average  does  not  strictly  represent  the  structure  in  any 
particular  location  because  of  the  lateral  variation  of  oce^ulic 
structure  with  the  age  of  the  sea  floor  (Forsyth,  1975). 


We  have  chosen  one  event  among  the  three  for  which  the 
data*  sasple  paths  predominantly  within  the  relatively  old 
(80-90  xn.y.  average  age)  central  Pacific.  The  locations  of 
the  event  (J^ril  26,  1973,  20*'26®30. 8® ,  latitude  19.9“N, 
longitude  155. 13 "W,  »  6.0),  stations  and  paths  eu:e  shown 

in  Figure  1. 

The  corresponding  phase  velocity  curve  is  taken  to  be 
a  weighted  average  of  'pure  path'  velocities,  calculated  using 
the  published  magnetic  anomalies  (Pitman  ^  al.  1974)  for  the 
Pacific  and  the  results  of  Forsyth  (1975,  1977)  on  the 
veuriation  of  phase  velocity  with  sea  floor  age.  We  divided 
the  Pacific  region  into  eleven  age  zones  plus  conti¬ 
nents  (Table  1) .  Eachr  great  circle  path  from  event  to 
station  was  plotted  on  the  map  of  age  zone  boundaries,  from 
%^ich  we  calculated  the  total  path  length  in  each  age  group 
for  the  sum  of  all  the  paths .  Then  the  weighted  average  of 
phase  velocity  and  group  velocity  at  each  period  were 
determined  using  the  pure-path  results  of  Forsyth  (1975,  1977), 
interpolating  where  necessary,  and  the  weights  shown  in  the 
Table  1.  The  resultauit  weighted  phase  velocity  c\irve  for 
the  central  Pacific,  appropriate  to  the  data  set  of 
Mitchell  et  al.  (1976) ,  Is  shown  in  Figure  2. 

The  magnetic  amomaly  map  on  the  basis  of  which 
seafloor  ages  were  estimated  does  not  include  the  marginal 
basins  of  the  western  Pacific,  across  which  pass  numy  of 


30. 


the  surface  wave  paths  used  here.  Additional  age  information 
for  these  basins  have  been  taken  from  Weissel  (1977)  for 
the  Lau  Basin,  from  Weissel  ^  al.  (1977)  for  the  Coral  Sea 
and  New  Hebrides  Basin,  from  Watts  et  al.  (1977)  for  the 
south  Fiji  Basin  and  from  Sclater  et  al.  (1976)  for  the 
Philippine  Sea. 

INVERSIONS 

Love  and  Rayleigh  waves  in  western  North  America.  The  starting 

models  for  calculation  of  partial  derivatives  for  western 

North  America  are  listed  in  Table  2  for  Q  independent  of 

frequency.  The  models  for  6  and  are  taken  from  Peqper  I. 

For  Q  independent  of  frequency  (and  for  other  dispersion- 

attenuation  relations)  the  starting  models  are  determined  by 

averaging  the  extreme  models  of  9  emd  eis  functions  of  depth. 

The  density  model  remains  real  but  is  allowed  to  vary  in  the 

Inversion  process.  The  P-wave  velocity  (Tedile  2)  is  chosen 

to  have  features  similar  to  model  CIT  111  of  Archambeau  et  al. 

(1969)'  euid  is  held  fixed  during  the  Inversion,  since  surface 

wave  phase  velocity  is  least  sensitive  to  changes  in  a  eunong 

(a,3,p)  and  thus  the  resolution  of  the  Inverse  problem  is 

poor  for  a.  Because  of  the  large  uncertainties  in  measured 

Q  we  will  assume  that  losses  under  purely  compressive 

—1  2  —1 

stress  are  negligible,  i.e.  ■  4/3(B/o)  Qg  .  Finally, 
a  2  km  sedimentary  layer  has  been  added  to  in^rove  the  fit 
for  Rayleigh  wave  phase  velocity.  For  purposes  of  calculating 
partial  derivatives,  the  equivalent  flat  model  to  that  in 
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Table  2  was  divided  into  28  layers  and  an  underlying  half space. 

The  cross-partial  elements  of  the  partial  derivative 

p 

matrix  (see  equation  7)  for  Rayleigh  waves,  9c^/382»  an 
order  of  magnitude  smaller  than  those  £or  Love  waves 
(Figure  3;  c.f..  Paper  I).  This  result  means  that  the 
difference  between  simultaneous  inversion  and  the  correction 
technique  of  Anderson  et  al.  (1977)  and  others  is  less 
signific2tnt  for  Rayleigh  waves  than  for  Love  waves.  However, 
the  advantage  of  simultaneous  inversion  still  remains  because 
the  chemges  in  the  model  itself  are  sensitive  to  the 
result  of  the  velocity  model  inversion. 

Resolution  analysis  was  conducted  as  described  above.  The 
optimal  averaging  kernels  (equation  15)  for  both  Love 
and  Rayleigh  waves  are  shown  in  Figure  4  for  five  layers. 

As  in  Paper  I,  a  5-layer  model  la  adopted  for  the  inversions. 

External  inversion  has  been  c^u^:ied  out  using  the  linear 
programming  technique  adopted  in  Lee  and  Solomon  (1975)  and  in 
Paper  t.  Envelopes  in  model  space  of  8^,  $2*  P  deter¬ 
mined  from  the  observational  data  and  their  uncertainties  aad 
from  the  constraints  that  the  n»del  parameters  lie  within 
certain  a  priori  bounds.  During  inversions,  the  shear  velocity 
and  attenuation  and  the  density  in  the  top  sedimentary  layer 
and  in  the  half  space  have  been  held  fixed. 

The  envelopes  of  shear  velocity  8,  attenuation  ^nd 

density  p  are  given  in  Table  3  and  Figure  5  for  a  model  (S31W) 
with  Q  independent  of  frequency,  and  for  models  (S32W,  S33W)  with 
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Q  vairying  as  a  power  of  the  frequency  (v  «  1/5  and  1/2,  respec¬ 
tively)  .  The  corresponding  envelopes  in  data  space  are 
illustrated  for  these  models  in  Figures  6-9.  A  1  Hz  reference 
frequency  is  used  for  all  models. 

Love  and  Rayleigh  waves  in  east-central  North  America. 

The  starting  models  of  <x,6,  P  and  for  east-central  North 

America  are  listed  in  Taible  4  .  The  models  of  a,  0  emd  P  are 

adapted  from  the  results  of  McEvilly  (1964)  and  the  model 

for  the  cons tan t-Q  Inversion  is  from  Lee  and  Solomon  (1975) . 

The  Initial  model  for  inversion  using  a  power  law 

dependence  of  Q  on  frequency  is  chosen  to  be  an  acceptable 

solution  to  the  separate  Q  inversion. 

The  averaging  kernels  Ej^j  from  resolution  analysis  for 

both  Love  and  Rayleigh  waves  are  shown  in  Figure  10,  for  five 

layers.  Four  layers  and  a  fixed-parameter  half space  were  used 

0 

in  the  inversions. 

The  envelopes  of  sheeu:  velocity  B,  attenuation  Qg^  and 
density  P  are  given  in  Table  5  and  Figure  11  for  a  dispersion 
model  (S31E)  with  Q  Independent  of  frequency  emd  for  a  model 
(S32E)v  with  Q  varying  as  a  power  of  frequency  (vb1/2)  .  The 
corresponding  envelopes  in  data  space  are  ill\istrated  for  these 
models  in  Figures  12-15. 

Rayleigh  waves  in  the  central  Pacific.  The  starting  models  of 
a,8>  P  and  Qg^  are  listed  in  Table  6.  The  starting  model  for 
the  crust  is  modified  from  the  'standard  crustal  section'  of 
Forsyth  (1975),  in  which  5  km  of  water  layer,  0.2  km  of 
sedimentary  layer  and  6.8  km  of  crustal  layer  are  assumed.  For 
the  starting  model  of  the  mantle,  the  density  is  taken  to  be 
3. 4-3. 5  q/ca? ,  the  S-wave  velocity  6  is  4.35  lOQ/sec  in  the  low 
velocity  zone  (50-220  km)  and  4.60  ka/sec  in  the  high  velocity  lid 
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a  'V'  1.7B  (<220  km  depth)  and  a  1.83  (>220  km  depth).  The 
starting  model  for  is  taken  from  Mitchell  (1976).  The 

p 

averaging  kernels  for  Rayleigh  waves  are  shown  in  Figure  16 
for  four  layers. 

The  envelopes  of  shear  velocity  3  and  sheeu:  attenuation 
are  given  in  TaOale  7  and  in  Figure  17  for  a  dispersion 
iBodel  (S21P)  with  Q  independent  of  frequency.  No  test  for  the 
frequency  dependence  of  Q  has  been  performed  since  the 
data  are  relatively  poor.  The  corresponding  envelopes  in  data 
space  for  model  S21P  aire  illustrated  in  Figures  2  and  18.  The 
dashed  lines  in  Figure  17  represent  the  average  of  the  models 
defining  the  envelopes  in  model  space.  This  averaged  model 
provides  a  good  fit  to  the  data  (Figures  2  and  18) . 

DISCUSSION 

Love  and  Rayleigh  waves  in  western  North  America.  A  problem  in 
western  North  America  may  be  SV/SH  emisotropy.  Love  and  Rayleigh 
wave  phase  velocities  are  incompatible,  in  the  sense  defined  by 
Lee  and  Solomon  (1975) ,  for  almost  the  entire  common  frequency 
range' (Figures  6  and  7) .  Also  as  noted  by  Lee  and  Solomon  (1975) , 
Rayleigh  wave  attenuation  in  the  period  range  35-40  seconds 
shows  disagreement  with  Love  wave  attenuation  This  is  not  a 

consequence  of  simultaneous  Inversion  but  appears  to  be  a 
consequence  either  of  auiisotropy  or  of  some  interference 
effects  in  the  waveforms  amalyzed.  The  measurements  of 
Rayleigh  wave  phase  velocity  in  Solomon  (1971)  aire 
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comparatively  lower  than  other  reported  meastirements  in  the 
western  United  States  (Biswas  and  Knopoff ,  1974) .  This  may  be 
a  real  effect,  however;  the  contour  map  of  velocity  in  the 
western  United  States  (Archaunbeau  et  ,  1969)  shows  that  the 
path  between  LON  2uid  TUG  lies  in  a  low  velocity  ' valley ' . 

In  Paper  I,  the  asstnnption  of  constemt  Q  led  to  the  removal 
of  the  requirement  for  a  low  velocity  zone  for  sheeu:  at  frequencies 
at  or  above  1  Hz.  When  Rayleigh  as  well  as  Love  wave  data  are 
used,  a  low  velocity  zone  at  1  Hz  is  a  feature  of  all  models 
resulting  from  the  inversion.  The  low  velocity  zone  and  low  Q 
zones  coincide,  a  property  that  is  preserved  even  when  st2u:ting 
models  with  no  low  velocity  zones  are  used  in  the  Inversion 
(Paper  I) . 

Love  and  Rayleigh  waves  in  east-central  North  America.  The 
simultaneous  inversion  results  for  east-central  North  America 
give  a  thick  lithosphere,  or  high  velocity  and  high  Q  lid,  and 
a  monotonically  increasing  velocity  model  at  1  Hz.  A  low 
velocity  zone  may  not  be  required  at  1  Hz.  However,  a  low  Q 
zone  is  probably  present  at  depths  greater  than  130  km. 

For  both  sets  of  North  Amerlceui  attenuation  data,  much  of 
the  data  Inconpatibility  in  certain  period  remges  which  occurred 
in  the  separate  inversion  for  Q  ^  (Lee  and  Solomon,  1975)  does 
not  arise  in  the  slmultameous  inversion.  This  is  true  in 
particular  of  the  low  values  near  20  to  25  sec  period,  a 
feature  of  other  data  sets  (Tryggvason,  1965;  Tsai  emd  Aki, 
1969) . 

Another  result  of  simultaneous  inversion  for  both  regions 
is  that  the  S-wave  velocities  in  the  mantle  are  Increased 
considerably  and  density  is  decreased  compared  to  the  results 
of  separate  inversion  of  phase  velocity.  At  the  same  time. 
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for  Q  independent  of  frequency,  Qg  in  the  asthenosphere  is 
greater  than  predicted  by  the  results  of  separate  inversion. 

The  frequency  dependent  models  are  also  satisfactory 
for  predicting  observations.  To  discriminate  eunong  the 
possible  frequency  dependencies  of  more  precise 

measiirements  of  or  of  body  wave  dispersion  must  be  made. 
Rayleigh  waves  in  the  central  Pacific.  The  res\ilts  of 
simultauieous  inversion  for  6  emd  Qg^  in  the  central  Pacific 
are  shown  in  Figxire  17.  The  low  velocity  zone  2uid  low  Q  zone 
coincide  and  are  extensive  (60  to  225  km  depth) .  The 
lithospheric  thickness  is  less  than  in  both  western  North 
America  (tectonic)  and  east-central  North  America  (stable) . 

The  result  of  simultaneous  inversion  (model  S21P)  shows  a 
shear  velocity  profile  generally  similau:  to  Forsyth's  (1977) 
old  ocean.  While  Forsyth  suggested  70-90  km  as  the  starting 
depth  of  partial  melting,  model  S21F  gives  a*  slightly  shallower 
asthenosphere  at  60  km  depth,  in  agreement  with  Mitchell  (1976) . 
This  may  be  due  to  the  influence  of  Q  ^  on  the  velocity  structure 
but  is  at  least  as  likely  due  to  the  coarse  layering  euid 
age-averaged  data  iised  in  this  study. 

Compared  with  Mitchell's  (1976)  Q  model,  the  noticeable 
differences  in  S21P  are  that  the  low  Q  zone  may  extend  deeper 
than  220  km,  though  the  Q  structure  at  that  depth  is  not  well 
resolved,  and  that  the  low  Q  zone  is  more  pronounced.  At 
around  100  sec  period,  Mitchell's  data  were  incompatible. 
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CONCLUSIONS 

We  have  extended  the  formalism  (Paper  I)  for  the 
simriltemeous  inversion  of  surface  wave  phase  velocity  and 
attenuation  to  Rayleigh  waves.  Simultsuieous  inversions  of 
data  for  both  Love  and  Rayleigh  waves  have  been  accomplished 
for  tectonically  active  aind  stable  continental  paths  euid  an 
inversion  for  Rayleigh  wave  data  alone  has  been  given  for  a 
at^d>le  ocean  basin. 

An  inportant  aspect  of  the  simultameous  inversion  is  the 
incorporation  of  the  intrinsic  auielastic  dispersion  which 
results  from  linearity  and  causality.  This  intrinsic 
dispersion  of  shear  velocity  can  be  quite  large  within  the 
low-Q  zone  beneath  oceanic  and  tectonically  active  regions. 
Beneath  western  North  America,  for  instance,  dispersion  can 
be  as  great  as  10  percent  over  two  decades  in  frequency, 
depending  on  the  frequency  dependence  of  Q 

The  simultaneous  Inversion  technique  presented  here  and 
in  Paper  I  is  formally  different  from  the  approximate  inversion 
scheme  of  Anderson  et  al.  (1977)  and  gives  different  results. 
These  differences  can  be  significant  when  the  observational 
frequency  band  is  wide  and/or  the  reference  frequency  is  far 
from  the  observational  frequencies.  No  matter  how  the 
observational  frequency  band  is  spread,  no  matter  where  we 
choose  a  reference  frequency,  a  given  dispersion- attenuation 
relationship  may  be  used  to  extrapolate  models  to  any 
frequencies.  The  difference  between  the  simultaneous 


inversion  approach  and  the  simple  correction  of  phase  velocity 
for  anelastic  dispersion  is  larger  for  Love  waves  them  for 
Rayleigh  waves,  because  the  cross  partial  derivatives  for 
Rayleigh  waves  are  much  smaller  than  for  Love  waves.  However, 
most  of  the  advemtages  of  simulteuieous  inversion  are  still 
retained  for  Rayleigh  waves. 

Other  consequences  of  simultaneous  inversion  are  that 
the  incompatibility  in  the  attenuation  data  set  decreases 
substemtially  emd  that  the  resolution  in  versus  depth  in 
the  earth  is  inproved  for  a  given  dispersion-attenuation 
relation  over  the  sepeurate  inversion  of  Q~^  data  alone. 

Several  possible  Intrinsic  dispersion-attenuation 
relations  were  tested  against  the  Love  and  Rayleigh  wave 
phase  velocity  and  data  sets  for  North  American  paths, 
including  those  for  Q  independent  of  frequency  ‘and  Q  varying 
as  a  small  positive  power  of  frequency.  The  availeible  surface 
wave  data  sets  do  not  discriminate  among  these  relations. 
Probably  the  best  prospect  for  choosing  aunong  them  is  in 
observing  body  wave  dispersion  (Jordu  and  Sipkln,  1977)  or 
frequency  dependence  of  Q  (Solomon  et  al. ,  1970;  Der  emd 
McElfresh,  1977). 

The  most  interesting  results  of  the  inversions  are  the 
lateral  variations  of  mantle  structure.  A  distinctive  low  Q 
zone  appears  .to  exist  everywhere  we  have  data,  including 
western  North  America  (tectonically  active  continent) ,  east- 
central  North  America  (stable  platform)  emd'the  central  Pacific 
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(stadsle  oceem  basin) .  However,  the  thickness  of  the  high-Q 
lid  varies  from  place  to  place:  60+20  km  in  the  central 
Pacific,  80+20  km  in  western  North  America  and  130+30  km  in 
east-central  North  America.  These  results  are  related  to  the 
differing  tectonic  history  of  each  region.  In  east-central 
North  America,  a  thick  lithosphere  has  presumably  grown  over 
time,  and  the  asthenosphere  shows  relatively  mild  attenuation 
and  admits  the  possibility  of  no  low  velocity  zone  at  a 
frequency  of  1  Hz  and  above.  In  western  North  America,  a 
thinner  lithosphere  with  substantially  higher  attenuation  in 
the  asthenosphere  is  characteristic.  Many  authors  have 
suspected  there  is  svibstantial  pairtial  melting  in  the 
asthenosphere  of  this  region.  The  cuccess  heat  necessary  to 
produce  melting  may  have  been  related  to  the  recent  subduction 
of  oceanic  lithosphere  along  western  North  America. 

The  low  velocxty  2md  low  Q  zones  coincide  in  western 
North  ^erica.  The  inversion  of  Love  wave  data  alone  (Paper  I) 
shows  that  the  assuoption  of  Q  independent  of  frequency  over 
the  coitire  seismic  band  leads  to  the  removal  cf  the 
requirement  for  a  low  velocity  zono  for  she^lr  waves  at 
freqxiencies  ed)Ove  1  Hz.  Ihe  low  velocity  zone  persists  at 
these  frequencies,  however,  in  the  results  of  the  combined 
inversion  of  Love  and  Rayleigh  wave  data.  A  possible 
anisotropy  in  western  North  America  may  have  contributed  to 
this  result,  though,  as  well  as  to  em  inability  to  discriminate 
among  various  dispersion-attenuation  relations.  In  the  Pacific 
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the  low  velocity  2uid  low  Q  zones  coincide  as  in  western  North 
America.  The  Pacific  shows  a  thinner  lithosphere  (  60  km  thick) 
but  lower  shear  attenuation  Q  ^  in  the  asthenosphere  than  in 

p 

the  tectonically  active  western  North  America. 
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APPENDIX 

PARTIAL  DERIVATIVES  AT  A  REFERENCE  FREQUENCY  FOR  RAYLEIGH  WAVES 

Using  equation  (28),  partial  derivatives  at  a  reference 
frequency  can  be  given  as  follows.  For  convenience,  the 
symbols  B°,  a°,  emd  a®  ®2'  “l'  ®2' 

respectively,  at  the  reference  frequency. 
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The  last  steps  of  equation  (A.  3)/  (A. 4),  (A.7)  »  'and  (A. 8) 
are  consequences  of  the  Cauchy~Rleinann  relations  (11). 
Front  dispersion  relations  (21)  emd  (26) 

2 

n  ln((ij/w  ) 
o  . 

gi  » 

irv  -V 

cot  —  tl-(<»j/Uo)  ] 


for  constant  Q 

for  power  law  Q 


for  constant  Q 


for  power  law  Q 


With  the  inversion  at  a  reference  frequency >  the  matrix 
equation  (6)  for  Love  waves  is  as  given  in  Paper  I.  For 
Rayleigh  waves,  the  matrix  equation  (7)  becomes 


where 


and 
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Table  1.  Fraction  of  great  circle  paths  in  each 
age  group  for  the  26  April  1973  Pacific 
event 


Path 

Distance  (*) 

Percentage 

0-5  m.y.  ocean 

19.99 

1.54 

5-10  " 

33.96 

2.61 

10-20  " 

98.86 

7.59 

20-38  " 

26.06 

2.00 

38-53  " 

205.87 

15.81 

53-65  " 

137.22 

10.54 

65-83  " 

23.00 

1.77 

83-100  • 

303.73 

22.33 

100-135  " 

179.77 

13.81 

135-190  " 

150.17 

11.53 

>  190  " 

5.86 

0.45 

North  Americam* 
continent 

32.57 

8.27 

South  American 
continent 

9.83 

0.76 

Total 

1302.04 

100 

*Continental  paths  in  southeast  Asia  are  included 
in  this  group. 
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Table  2.  Starting  model  for  inversion  of  combined  Love 
and  Rayleigh  phase  velocity  wd  attenuation  in 
western  North  America  (case  v  »  o) 


Depth,  km 

8 ,  km/sec 

a,  km/sec 

p,  g/cm^ 

100/Q 

0-2* 

1.72 

4.00 

2.21 

.32 

2-13 

3.58 

6.00 

2.80 

.32 

13-21 

3.58 

6.20 

2.80 

.32 

21-45 

3.96 

6.50 

3.25 

.13 

45-64 

4.54 

7.619 

3.40 

.13 

00 

1 

4.54 

7.615 

3.40 

.13 

84-94 

4.48 

7.615 

3.45 

5.93 

94-128 

4.49 

7.619 

3.45 

5.93 

128-160 

4.49  . 

7.622 

3.45 

5.93 

160-180 

5.27 

7.90 

3.50 

2.96 

180-220 

5.27 

8.23 

3.50 

2.96 

220-260 

5.27 

8.25 

3.50 

2.96 

260-300 

5.27 

8.27 

3.50 

2.96 

300-350 

5.27 

8.43 

3.50 

2.96 

Halfspace* 

5.27 

8.53 

3.90 

2.96 

*Paraineters  fixed  during  inversion  for  these  layers. 
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Table  3.  Envelopes  of  shear  velocity  and  sheeur  attenuation 
at  1  Hz  euid  of  density  front  the  simultaneous 
inversion  of  Love  euid  Rayleigh  wave  phase  velocity 
and  attenuation  in  western  North  America. 


Layer 

no. 

Depth 

(km) 

minimum  B,  km/sec 

S31W  S32H  S33W 

maximum 

S31W 

8  /  km/sec 

S32W  S33W 

2 

2-21 

3.527 

3.529 

3.533 

3.492 

3.621 

3.637 

3 

21-45 

3.897 

3.886 

3.812 

4.012 

4.052 

4.029 

4 

45-84 

4.382 

4.259 

4.293 

4.583 

4.549 

4.690 

5 

84-160 

4.105 

4.012 

3.823 

4.463 

4.412 

4.263 

6 

160-350 

4.782 

4.750 

4.702 

5.500* 

5.500* 

5.500* 

Layer 

no. 

Depth 

(km) 

minimum  100/Qg 

S31W  S32W  S33W 

maximum 

S31W 

100/Qfl 

S32H 

S33H 

2 

2-21 

0.274 

0.148 

0.051 

0.490 

0.282 

0.118 

3 

21-45 

0.059 

•0.0 

0.0 

0.250 

0.149 

0.105 

4 

45-84 

0.0 

0.0 

0.0 

0.447 

0.137 

0.087 

5 

84-160 

4.200 

2.000 

0.453 

5.330 

2.337 

0.783 

6 

160-350 

1.350 

0.119 

0.084 

3.640 

0.782 

0.860 

Layer 

no. 

'  Depth 
(km) 

minimum  p ,  g/cm^ 

S31W  S32W  S33W 

maximum 

S31W 

p,  g/cm 
S32W 

3 

S33W 

2 

2-21 

2.807 

2.846 

2.850 

2.90* 

2.90* 

2.90* 

3 

21-45 

3.070 

3.042 

3.000* 

3.30* 

3.30* 

3.30* 

4 

45-84 

3.200* 

3.200* 

3.200* 

3.434 

3.50* 

3.50* 

5 

84-160 

3.200* 

3.200* 

3.200* 

3.491 

3.441 

3.50* 

*A  priori  bounds  not  well  resolved  by  the  inversion 
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Table  4.  Starting  model  for  inversion  of  Love  emd 

Rayleigh  wave  phase  velocity  euid  attenuation 
in  east-central  North  America  (case  v  =  0) 


Depth,  km 

8,  km/sec 

a,  km/sec  . 

p,  g/aa? 

100/Q| 

0-11 

3.50 

6.10 

2.9 

0.1 

11-20 

3.68 

6.20 

2.9 

0.1 

20-38 

3.94 

6.40 

2.9 

0.1 

38-62 

4.75 

8.15 

3.3 

0.1 

62-102 

4.61 

8.20 

3.3 

0.1 

102-135 

4.45 

8.20 

3.4 

0.1 

135-350 

4.45 

8.20 

3.4 

2.3 

Halfspace* 

4.80 

8.70 

3.6 

2.3 

*Peu:euneters  fixed  during  inversion  for  this  layer 


Table  5.  Envelopes  of  shear  velocity  and  shear  attenuation 

at  1  Hz  and  density  from  the  simultaneous  inversion 
of  Love  and  Rayleigh  wave  phase  velocity  emd 
attenuation  in  east-central  North  America 


Layer 

no. 

Depth 

(km) 

Minimum  B 
S31E 

,  Icm/sec 
S32E 

Maximum  B  > 
S31E 

r  km/sec 
S32E 

1 

0-38 

3.684 

3.659 

3.692 

3.666 

2 

38-82 

4.557 

4.478 

4.594 

4.513 

3 

82-134 

5.092 

4.719 

5.321 

4.890 

4 

134-350 

5.200 

5.147 

5.500* 

5.500* 

Layer 

no. 

Depth 

(km) 

Minimum  100/Qg 

S31E  S32E 

Maximum  lOO/Qg 
S31E  S32E 

1 

0-38 

0.165 

0.037 

0.263 

0.055 

2 

38-82 

0.0 

0.0 

0.222 

0.048 

3 

82-134 

3.056 

0.241 

5.630 

0.831 

4 

134-350 

0.182 

0.036 

3.846 

0.364 

Layer 

no. 

Depth 

(km) 

Minimum  p 
S31E 

,  g/cm^, 
S32E 

Haximian  P  r 
S31£ 

g/cm^ 

S32E 

1 

0-38 

2.896 

2.895 

2.900* 

2.900* 

2' 

38-82 

3.200* 

3.000* 

3.241 

3.010 

3 

82-134 

3.200* 

3.200* 

3.327 

3.220 

4 

134-350 

3.600* 

3.600* 

3.800* 

3.800* 

*A  priori  bounds  not  well  resolved  by  the  Inversion. 
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Table  6.  Steurting  model  for  the  inversion  of  Rayleigh  wave 
phase  velocity  emd  attenuation  in  central  Pacific 


Depth f  km  Sr 

km/sec 

a ,  km/sec 

p,  q/aa} 

100/Qg 

0-5* 

0.0 

1.52 

1.03 

0.0 

5-5.2* 

1.0 

1.70 

2.20 

0.80 

5.2-12* 

3.76 

6.53 

3.00 

0.60 

12-47 

4.60 

8.17 

3.34 

0.60 

47-67 

4.35 

7.70 

3.40 

0.50 

67-107 

4.35 

7.70 

3.44 

1.00 

107-147 

4.35 

7.70 

3.44 

1.20 

147-227 

4.35 

7.70 

3.44 

1.00 

227-302 

4.60 

8.49 

3.50 

0.60 

302-402 

0 

4.80 

8.81 

3.50 

0.40 

Half space* 

5.00 

9.00 

3.76 

0.40 

*Paraiiieters  held  fixed  during  inversion  for  these  layers. 
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Table  7.  Envelopes  of  shear  velocity  and  shear  attenuation 
at  1  Hz  in  the  Pacific 


Layer  no. 

Depth  (loa) 

Minimum  8 

Maximum  6 

4 

12-57 

4.542 

4.571 

S 

57-107 

4.195 

4.350 

6 

107-227 

4.212 

4.516 

7 

227-402 

4.50* 

5.00* 

Layer  no. 

Depth  (Ion) 

Minimum  100/Qg  Maximum  lOO/Qg 

4 

12-57 

0.413 

0.440 

5 

57-107 

1.295 

1.742 

6 

107-227 

0,900 

1.695 

7 

m 

227-402 

0.404 

2.000 

*A  priori  bounds  not  well  resolved  by  the  inversion. 
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Figure  Captions 

Fig..  1.  Pacific  area  map  showing  paths  between  the 

i^ril  26,  1973  earthquake  and  stations  (ALQ,  ANP, 

ARE,  BAG,  BOG,  CHG,  COL,  COR,  DAV,  HKC,  JCT,  LPB, 

NIL,  RAR,  RIV,  SNG,  TAU,  TUC,  WEL)  used  in  the  q"^ 
measurements  of  Mitchell  et  al.  (1976).  Azimuthal 
equidistauit  projection  about  0**N,  180°E. 

Fig.  2.  Rayleigh  wave  phase  velocity  for  the  central 

Pacific.  Circles  with  error  baurs  are  calculated  as  a 
weighted  average  of  'pure  path*  velocities  (Forsyth, 

1975,  1977) ,  including  measurement  uncertainties,  using 
the  magnetic  amomaly  map  of  Pitman  et  al.  (1974).  The 
solid  lines  are  the  envelope  associated  with  extremal 
earth  models  for  Inversion  S21P.  The  dashed  line  is 
the  predicted  dispersion  for  the  model  shown  as 
dashed  in  Figure  17. 

Fig.  3.  Selected  partial  derivatives  of  the  real  part  of 
Rayleigh  wave  phase  velocity  with  respect  to  the 

p 

imaginary  part  of  shear  velocity  (9c£/96  2'  solid  lines)  and 
to  the  real  part  of  shear  wave  velocity 

Ocf/36]^,  dashed  lines)  per  unit  layer  thickness  for 
the  initial  model  in  Table  2.  The  partials  shown  are 
for  frequency-independent  at  the  frequency 
Indicated;  for  frequency-dependent  Q  ^  relations,  the 
partial  derivatives  have  a  similar  structure. 
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Discontinuities  in  the  partials  occur  at  discontinuities 
in  the  initial  model. 

Fig.  4.  Resolving  kernels  for  (a)  shear  velocity  and  (b) 
shear  attenuation  at  selected  depths  (arrows)  at  the 
reference  frequency  1  Hz,  using  both  Love  and  Rayleigh 
wave  data  in  western  North  America  for  the  case  of  Q 
independent  of  frequency.  Model  standard  deviations  are 
shown  at  the  right. 

Fig.  5.  Envelopes  of  shear  velocity  emd  shear  attenuation, 
at  a  frequency  of  1  Hz,  and  density  for  SKsdels  S31W, 

S32W  and  S33W.  Solid  lines  represent  envelopes  of  S31W 
(v  a  0) ,  short-dashed  lines  are  for  S32W  (v  *  1/5)  and 
long-dashed  lines  eure  for  S33H  (v  *  1/2) . 

Fig.  6.  Love  wave  phase  velocity,  western  North  America. 
Observations  are  shown  by  circles;  vertical  bars 
represent  standard  deviations.  The  envelope  (solid 
l^Jies)  is  associated  with  the  extreme  earth  model 
bounds  from  inversion  S31W.  Open  circles  are  Incompatible 
data  for  this  inversion. 

Fig.  7.  Rayleigh  wave  phase  velocity,  western  North  America. 
Observations  are  shown  by  circles;  vertical  baurs 
represent  standard  deviations.  The  envelope  (solid  lines) 
is  associated  with  the  extreme  earth  model  bounds  from 
inversion  S31W.  Open  circles  are  incompatible  data  for 
this  inversion. 
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Fig.  8.  Love  wave  attenuation*  western  North  America. 

Circles  are  observations;  vertical  beurs  represent 
steuidard  deviations.  The  solid  lines  represent  the 
envelope  associated  with  extremal  earth  models  for 
inversion  S31W. 

Fig.  9.  Rayleigh  wave  attenuation*  western  North  America. 
Circles  are  observations;  vertical  bars  represent 
standard  deviations.  The  solid  lines  represent  the 
envelope  associated  with  extremal  earth  models  for 
inversion  S31W.  Open  circles  are  incoaqpatible  data  for 
this  inversion. 

Fig.  10.  Resolving  kernels  for  (a)  shctar  velocity  and  (b) 
shear  attenuation  at  selected  depths  (aurrows)  at  the 
reference  frequency  1  Hz,  using  both  Love  and  Rayleigh 
wave  phase  velocity  amd  attenuation  in  east-central  North 
America  for  the  case  of  Q  independent  of  frequency. 

Model  standaurd  deviations  aure  shown  at  the  right. 

Fig.  11.  Envelopes  of  sheau:  velocity  emd  sheau:  attenuation, 
at  a  frequency  of  1  Hz,  and  density  for  models  S31E  and 
S32E.  Solid  lines  represent  envelopes  of  S31E  (v  <■  0), 
auid  dashed  lines  aure  for  S32E  (v  «  1/2). 

Fig.  12.  Love  wave  phase  velocity,  east-central  North  America. 
Observations  are  shown  by  circles;  vertical  bars  represent 
standard  deviations.  The  envelope  is  associated  with  the 
extreme  earth  model  bounds  from  inversion  S3 IE.  Open 
circles  are  incompatible  data  for  this  inversion. 


Fig.  13.  Rayleigh  wave  phase  velocity,  east-central  North 
America.  Observations  are  shown  by  circles;  vertical 
bars  represent  stemdard  deviations.  The  envelope  is 
associated  with  the  extreme  earth  model  bounds  from 
Inversion  S31£.  Open  circles  are  incompatible  data  for 
this  inversion. 

Fig.  14.  Love  wave  attenuation,  east-centred.  North  America. 

Circles  are  observations;  vertical  bars  represent  standard 
deviations.  The  solid  lines  represent  the  envelopes 
associated  with  extremal  earth  models  for  inversion  S31E. 

Open  circles  are  incompatible  data  for  this  inversion. 

Fig.  15.  Rayleigh  wave  attenuation,  east-central  North  America. 
Circles  are  obsezfvations;  vertical  bars  represent  standard 
deviations.  The  solid  lines  represent  the  envelope 
associated  with  extremal  euth  models  for  inversion  S31E. 

Fig.  16.  Resolving  kernels  for  (a)  sheeir  velocity  and  (b)  shear 
attenuation  at  selected  depths  (arrows) ,  using  Rayleigh 
wave  data  in  the  central  Pacific  for  the  case  of  Q 
Independent  of  frequency.  Model  standard  deviations  are 
shown  at  the  right. 

Fig.  17.  Envelopes  of  she^u:  velocity  (at  1  Hz)  and  shear 
attenuation  for  the  central  Pacific,  model  S21P  (Q~^ 
independent  of  frequency) .  The  dashed  lines  show  the  average 
of  the  models  making  up  the  envelopes. 


Fig.  18.  Rayleigh  wave  attenuation,  central  Pacific.  Circles 
are  observations;  vertical  bars  represent  stemdard 
deviations.  The  solid  lines  represent  the  envelope 
associated  with  extremal  earth  models  for  Inversion  S21P. 
Open  circles  are  incompatible  data  for  this  Inversion. 

The  dashed  line  is  predicted  for  the  model  shown  as 
dashed  in  Figure  17. 
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